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We discuss the phase structure of dense matter, in particular the nature of the transition between 
hadronic and quark matter. Calculations within a Ginzburg-Landau approach show that the axial 
anomaly can induce a critical point in this transition region. This is possible because in three-flavor 
quark matter with instanton effects a chiral condensate can be added to the color-flavor locked 
(CFL) phase without changing the symmetries of the ground state. In (massless) two-flavor quark 
matter such a critical point is not possible since the corresponding color superconductor (2SC) does 
not break chiral symmetry. We study the effects of a nonzero but finite strange quark mass which 
interpolates between these two cases. Since at ultra-high density the first reaction of CFL to a 
nonzero strange quark mass is to develop a kaon condensate, we extend previous Ginzburg-Landau 
studies by including such a condensate. We discuss the fate of the critical point systematically and 
show that the continuity between hadronic and quark matter can be disrupted by the onset of a 
kaon condensate. Moreover, we identify the mass terms in the Ginzburg-Landau potential which 
are needed for the 2SC phase to occur in the phase diagram. 

PACS numbers: 12.38. Mh,24.85.+p 



I. INTRODUCTION 



Cold, strongly interacting matter appears as ordinary nuclear matter at small baryon chemical potential and 
as color-flavor locked (CFL) quark matter P, Q at sufficiently large baryon chemical potential. The former is an 
observational fact while the latter follows theoretically from first-principle calculations within perturbative QCD. 
Since QCD is asymptotically free [!, H| these calculations are valid at asymptotically large chemical potential. How 
matter evolves from nuclear matter to CFL quark matter, i.e., which phases and phase transitions it passes in between, 
is an unsolved problem. This problem is of theoretical relevance since it addresses a poorly understood territory of 
the QCD phase diagram. It is also of phenomenological relevance since the interior of compact stars contains such 
"in-between" matter, i.e., matter neither well described by low-density nuclear physics nor by ultra- high-density 
perturbation theory (for an introduction to dense matter in compact stars see Ref. [5[). 

From the symmetries of ordinary nuclear (not hypernuclear) matter and the CFL phase it is clear that there must 
be at least one true phase transition between these two phases. However, since CFL breaks chiral symmetry it is 
very similar to hadronic matter and in principle allows for the intriguing possibility of a quark-hadron continuity 
Q. In other words, if hadronic matter is "prepared" by a series of appropriate phase transitions, including onsets 
of hyperonic matter and subsequent transitions to hyperon supcrfluids, it can become indistinguishable from CFL 
quark matter 0] ■ Whether these transitions happen, in particular whether they happen before deconfinement sets in, 
depends on the value of the constituent strange quark mass and the details of the baryonic interactions. 

We may also think of the quark-hadron transition by starting from the high-density limit. Upon going down in 
density, the effect of the strange quark mass becomes more important, imposing a stress on the symmetric pairing 
pattern of CFL. Systematic studies at high densities show that CFL reacts on this stress first by developing a kaon 
condensat e followed by an anisotropic phase with a kaon supercurrent [9i-fTTI| . and then possibly a crystalline 
phase Also pairing pat terns w hich do not break chiral symmetry may arise, such as the 2SC phase jl5| or 

spin-one color superconductors [h6l - [T8j . Only if such a phase does not appear before confinement sets in, a smooth 
quark-hadron crossover is possible. 

The problem of the quark-hadron continuity can thus be thought of in the following illustrative way. Start from 
low (high) density and go up (down) in density and determine all phases on the way, say at zero temperature and for 
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the moment ignoring the dcconfincmcnt (confinement) transition. Although we currently do not have the theoretical 
tools to determine all these phases in detail, this will in principle lead to two sequences of phases. One describes 
strongly interacting hadronic matter between ordinary nuclear matter at low density and some phase at very high 
density, the other describes quark matter between CFL at ultra-high density and some phase at lower density. Now 
we would like to put these two sequences together to obtain the true QCD phase diagram at small temperatures. We 
know that they fit smoothly if we glue them together at the point of the CFL phase in the quark matter part and the 
point of the hyperon supcrfluid phase described in Ref. Q in the hadronic part. 

In this paper we address this problem in a very sim plifi ed fashion. Our first simplification is the use of a Ginzburg- 
Landau free energy as in Refs. Qjl, [2(|, see also Refs. }2lM26j ; for recent discussions of the quark- hadron transition in 
a Nambu-Jona-Lasinio (NJL) model see Refs. [27i - [29| . Within the Ginzburg-Landau approach we will not be able to 
make firm predictions for QCD. The benefit is, however, that our study is model-independent and systematic since 
the free energy is based on the underlying symmetries and we can in principle determine and classify all phases in 
the parameter space of the Ginzburg-Landau coefficients. The order parameters we use are the chiral condensate and 
the diquark condensate. Our "nuclear matter" is thus mimicked simply by a nonvanishing chiral condensate, and we 
cannot account for the various phase transitions which potentially "prepare" hadronic matter for the quark-hadron 
crossover. The onset of a CFL diquark condensate is a very crude version of these transitions. After this onset 
has happened, the question is whether the chiral condensate can approach zero smoothly, i.e., without causing an 
additional phase transition. For the case of a vanishing strange quark mass this is possible, and the Ginzburg-Landau 
parameter space where this is realized has been determined in Refs. [l9l [20|. In particular, it can only be realized in 
the presence of a term that violates the U(1)a symmetry and thus is an effect of the axial anomaly (other effects of 
the anomaly in high-density QCD are discussed for instance in Refs. [30l. l3lj). Only then do the CFL phase without 
and the CFL phase with chiral condensate have the same residual symmetry. As a consequence, the £/(l)A-violating 
term - which couples chiral and diquark condensates - induces a critical point in the Ginzburg-Landau phase diagram 
where a first-order phase transition line ends. The Ginzburg-Landau study also shows that for infinite strange quark 
mass, i.e., in two- flavor quark matter, there is always a true phase transition from the phase where both chiral and 
diquark condensates are nonzero to the phase where only the diquark condensate is nonvanishing. 

We know that the real world is somewhere in between, i.e., the strange quark mass is neither negligibly small 
nor approximately infinite in the region of intermediate densities we are interested in. Therefore, in this paper we 
extend the previous Ginzburg-Landau studies by including a nonzero but finite strange quark mass. The first obvious 
consequence is the appearance of additional terms in the Ginzburg-Landau free energy. Additionally, we take into 
account meson condensation in CFL. Since CFL breaks chiral symmetry, there is an octet of Goldstone modes, just 
as in the hadronic phase. At high density an effective theory for these modes, based on the symmetries of CFL, can 
be derived Q, and it can be shown that the masses of the Goldstone modes are inversely ordered compared to the 
ordinary mesons in hadronic matter [32[ . The high-density effective theory suggests that the neutral kaons form a Bose 
condensate in the presence of a nonzero strange quark mass Q. The resulting phase has been termed CFL- if (see 
for instance Refs. |33U35l ] for properties of CFL- if computed from the high-density effective theory and Refs. j3a - [39| 
for calculations within an NJL model). It is thus natural to include kaon condensation into the Ginzburg-Landau 
calculation. We note, however, that at lower densities the meson masses may receive large instanton corrections [3lj . 
More formally speaking, at ultra-high densities, where instanton effects arc suppressed, the meson masses squared 
are quadratic in the quark masses since a linear term in the masses is forbidden by the U(1)a symmetry. When a 
linear, [/(l^-violating term is allowed at lower densities, the masses - and thus also their ordering - is subject to 
potentially large corrections and it is unclear which of the mesons form a Bose condensate. Therefore, our choice of 
a kaon condensate is inspired by the high-density results, but in principle also other meson condensates should be 
considered. 

Our paper is organized as follows. In Sec. [IT] we set up the Ginzburg-Landau potential, including mass terms and 
kaon condensate. Then we first evaluate the potential in the limit case of vanishing kaon condensate. This is done in 
Sec. IIII1 The purpose of this section is to explain the Ginzburg-Landau phase diagram in terms of the symmetries of 
the various phases and in the presence of a mass-induced linear term in the chiral potential. We build on the phase 
diagrams obtained in this section when we include the effect of a kaon condensate in Sec. IIVI In Sec. [V] we discuss the 
appearance of the 2SC phase due to mass corrections in the Ginzburg-Landau potential before we give our conclusions 
in Sec. EH 

II. GINZBURG-LANDAU POTENTIAL 

We are interested in a Ginzburg-Landau potential for the chiral and diquark condensates, including a neutral kaon 
condensate and corrections from the strange quark mass, based on Refs. [l9|, [2(|. The final result of this section is 
Eq. and the following pages are devoted to the derivation of this equation. 
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A. Symmetries and order parameters 

We shall consider a Ginzburg-Landau free energy of the form 

n = + n d + n 9d , (i) 

with a chiral part f2$ depending only on the chiral condensate $, a diquark part depending only on the diquark 
condensates and d R (which in turn depend on the kaon condensate) and an interaction part f2$d which cou- 
ples $ with dh and d R . In terms of left- and right-handed quark fields qL and q R we have (q Ri qLj) oc $ 3 i and 

{qlfiq^j) oc —e a ^ A CijB[d\ j ]g, (q Ri C 'q R j) oc e a P A £ijB[d R ]g, with flavor indices i,j,B, color indices a, (3, A, and the 
charge conjugation matrix C = «7 2 7 . Since we consider a three-flavor system, $, d^, and d R are 3x3 matrices. The 
quark fields transform under the symmetry group 

G = SU{3) C x SU{3) L x SU(3) R x U(1) b x U(l) A (2) 

as 

q L -> e^e-^Vhtfffi , to -> e fafl e faA V fl ^?fl , (3) 

where (Vl,Vr) € SU(3)l X SU(3) R is a chiral transformation, U £ SU(3) C is a color gauge transformation, 
giQ B g J7(l)s is a transformation associated with baryon number conservation, and e tolA £ U(\)a is an axial trans- 
formation. In terms of left- and right-handed U(l) transformations we have qL — > e laL qL, q R — > e taR q R , hence the 
vector and axial U(l) transformations in Eq. ([3]) follow from a R — {olr + ol)/2, ol a = {ot R — az,)/2. Eventually, 
our potential will not be invariant under the full group G. The chiral group SU(3)l x SU(3) r and the axial U(1)a 
become approximate symmetries after including the effects of a small strange quark mass and the QCD axial anomaly, 
respectively. 

The transformation properties of the order parameters under G are 

$^ e- 2iaA V L ^V^, (4) 

and 

d L -> e- 2laB e 2laA V L d L U T , d R e~ 2laB e - 2taA V R d R U T . (5) 
The mass terms are generated by the field 

( m u \ 

M = m d , (6) 
\ m s J 

which transforms under G in the same way as the chiral field $, 

M^e^FiMy]. (7) 

Although we shall write down the Ginzburg-Landau terms with general quark masses m u , md, m s , we shall later 
neglect the up and down quark masses and only keep the strange quark mass. 
Our ansatz for the order parameters is as follows. The chiral field is 

( On \ 

$ = a d . (8) 

V o o <j s J 

We shall derive the potential Q within this general ansatz, but later set for simplicity a u = a d = <J S . Different values 
for each quark flavor are more realistic in the presence of a strange quark mass and a kaon condensate. However, 
this would introduce additional independent parameters into our potential, making a systematic evaluation very 
complicated. Therefore, we shall use the symmetric case as a simplification. 
For the diquark condensate we use the ansatz 

10 

d R = d\ cosO/2) isinO/2) ) , (9) 
isin(0/2) cos(<£/2) 
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where 4> is the kaon condensate. For cf> = we recover the pure CFL order parameter d,L = dn = diag(d, d, d). A 
nonzero <p introduces a relative rotation between left- and right-handed diquarks. We have chosen a neutral kaon 
condensate which is the most likely possibility at high densities but, as explained in the introduction, at intermediate 
densities there could also be other meson condensates. 

As a comparison, note that the chiral field in the high-density effective theory of CFL [§| is 



(10) 



with the meson fields 8 a , the Gell-Mann matrices A a , and the (CFL version of the) pion decay constant / w . Ignoring 
all meson fields other than the neutral kaon, we set all 6 a 's to zero except for 0§ and 8j. Without loss of generality we 
can also set 9r = 0. Then, our condensate <f> is the vacuum expectation value of 6%/ and the chiral field is written 
in terms of the diquark condensates as 



dLd r R 
d? 




(11) 



While d^ and dp are gauge variant quantities, £ is gauge invariant. Therefore, our ansatz ([9]) is one of infinitely many 
choices for d^ and d R - all related by gauge transformations - which lead to Eq. (|11[) . The chiral field transforms 
under G as 

d L d\^e UaA V L d L d^Vl. (12) 

This is the same transformation property as the ordinary chiral field except for the transformations under U(1)a- 
This difference reflects the fact that in CFL the mesons are composed of four quarks, not two. 



B. Chiral potential 

We can now derive the explicit form of the potential f2. For the chiral part fi$ we collect all terms up to fourth 
combined order in M and $ with at most one power in the mass field M. The terms of 0(M°) which are invariant 
under G are 

Tr[$t$] = a l +a j + a ^ (l 3a ) 

Cn[#t$])2 = ((T 2 +ff 2 +(7 2 )2j (13b) 

Tr[($t4>) 2 ] = CT 4 +(r 4 + CT 4, (13c) 
If we do not require the potential to be invariant under U(1)a, we have the additional term pol ] 

det $ + h.c. oc e abc ei jk $ ai $ bj $ ck + h.c. = I2a u a d a s , (14) 

which, because of 

det$->e- 6iaA det$, (15) 

is only invariant under the discrete subgroup 2^(6) C U(1)a, as expected from the anomaly. Microscopically, det$ 
accounts for an effective six-point instanton vertex which converts three left-handed quarks into three right-handed 
quarks, thus violating axial charge conservation by an amount 2Nf = 6. 

The terms of order 0(M 1 ) arise from replacing one chiral field $ by the mass field M in each of the above O(M ) 
terms. We obtain 

Tr[M$ f ] + h.c. = 2(m„er„ + m d a d + m s a s ) , (16a) 
T*[M$ t ] , R[# t *]+h.c. = 2(m u <T u + m d <r d + m s <j s )(al + a 2 d + a 2 s ) , (16b) 
Tr[Af$$ t $] + h.c. = 2{m u al + m d al + m s a 3 s ), (16c) 

and the anomalous term becomes 



£abc£ijkM ai $ b: j$ ck + h.c. = 4(a u a d m s + <r u a s m d + cr d a s m u ) . 



(17) 
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We can now add the contributions (|13|) . (fT4|) . (fT6|) . p7|) . each coming with a separate prefactor, to obtain the chiral 
potential. Approximating m u ~ m c i ~ and setting for simplicity a u = o& = cr s = a, we can write the potential as 

a 1 +m s a 2 2 . ci+m s c 2 3 6 4 , . 

= flops' 7 H ^ c H g 4°" ' ' ' 

In the given approximation, m s gives rise to a linear term in a and yields corrections to the quadratic and cubic 
terms. Because of the linear term, the chiral condensate cannot vanish exactly in the ground state. Instead of a 
vacuum phase with a = there will be a phase with very small a where chiral symmetry is approximately restored 
and which is continuously connected to the chirally broken phase in which a has a sizable value. This is the most 
obvious consequence of the mass term. Due to the coupling of chiral and diquark condensates, to be discussed in Sec. 
Ill Dt we shall find other, less obvious, effects of the linear term for our phase diagram. These effects are discussed in 

secnni 



C. Diquark potential and comparison to high-density effective theory 

For the diquark potential fid we also start from the terms up to 0(d i ), first without mass insertions. Within our 
CFL-A' ansatz they simply yield the structures d 2 and d 4 since the kaon condensate always drops out, 

Tr[d L d{] = Tr[d R d R ] = 3d 2 , (19a) 
(Tr[d L d{}) 2 = (Tr[d R d R }) 2 = Tr[d L d{}Tr[d R d R ] 
= 3Tr[(d L d[) 2 ] = 3Tr[(d R d R ) 2 } = 3Tr[d R d f L d L d R ] = ^ ■ ( 19b ) 

All these terms are invariant under the full group G. There is no such term as det dL yR since this term would not 
only break U{1)a but also baryon number conservation which must not be explicitly broken. To include the effect of 
quark masses, we need to replace di,d R by M. From Eqs. (fT9|) the only possible term (except for a term constant in 
dL, da which we can omit) is 

TT[d[Md B ] + h.c. = 2d 2 [?7i n + (m d + m s ) cos c/>] . (20) 

This term is invariant under Z^(6) C U(1)a and thus is allowed in the presence of the anomaly. Other 0(M 1 ) terms 
arise from replacing one chiral field dirf^ in the 0(d 6 ) terms. They all yield the same structure and are invariant 
under 2^(6), 

3(Tr[4MdR<4efo] + h.c.) = 3(Tv[d[Md R d[d L ] + h.c.) 

= Tr[4-/Wd R ]Tr[dJjd fl ] + h.c. = Tr[d[Md R }Tr[dld L } + h.c. = 6d 4 [m u + (m d + m s ) cos </>] . (21) 

So far the only structure we have produced for the kaon condensate is cos <j>. If we were to stop here, the minimization of 
fid would not allow for nontrivial condensates. This would be in contradiction to high-density calculations. Therefore 
we need to include at least one extra term with nontrivial structure in 0. To find this term it is useful to briefly 
discuss the Lagrangian of the high-density effective theory, 

£ cff = ^Tr[V £Vo£t - v 2 diY,ditf] + BTr[M£t + Aft£] + ^detMTrfM" 1 ^ + S f )] , (22) 
where weak-coupling calculations give the following values for the constants [1, [U HH , 

/i = ™^, „;=i, ._»{£, ...(i^an^i^. (23, 

Here, \x q is the quark chemical potential, A the superconducting gap parameter, c ~ 0.155, g the strong coupling 
constant, and Aqcd the QCD scale factor. The term linear in M with prefactor B corresponds to the term in Eq. (|2"0")) 
and breaks U{1)a- An important contribution to the effective Lagrangian is given by promoting the time derivatives 
to covariant derivatives, 

MAP 

V E = d Z + i[A,Y,], A = — . (24) 
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Including only neutral kaons, i.e., using the chiral field from Eq. (fTTj) . the treedevel potential V{4>) = — C e s becomes 
[after subtracting the "vacuum" contribution V((f> = 0)] 



V{4>) = - f -f^ 2 sin 2 + />^o(l - cos0) , (25) 



with the effective kaon mass squared 

2 _ 2B(md + m s ) . . 

m K o = — h am u (m s + rrid ) ■ (26) 

The kaon mass squared receives an instanton contribution linear in the quark masses [3ll . l4lj and a contribution 
quadratic in the quark masses. In the weak-coupling limit at very high densities, [i q ^ Aqcd, the instanton contri- 
bution is negligible and U(1)a is effectively restored. In this case the term M _1 detM becomes dominant, resulting 
in the inverse meson mass ordering (32j . We do not include a term of 0(M 2 ) in our Ginzburg-Landau potential, but 
nevertheless assume the kaon to be the most likely meson to condense. 

The covariant derivative (|24p gives rise to the term oc sin 2 cj) which contains the effective kaon chemical potential 

" s ^r- < 27 > 

This term is crucial for kaon condensation to occur. One can see this for instance by expanding V((f>) for small </>, 
say up to order </> 4 to reproduce the potential of an ordinary </> 4 model. The resulting <fi 2 term becomes negative only 
for fi > m^-o in which case a nonzero value of </> minimizes the potential. We thus conclude that the critical strange 
quark mass for the onset of kaon condensation scales at high density as m s ~ mi A 2 / 3 . 

From Eq. (fiM)) we see that the term oc /j, 2 in the potential is formally of order M 4 . Nevertheless, it can become 
comparable to the M 2 terms ||, and therefore we need to include this term in our potential, 1 

^[[MM^dLd^M^M.dRd^]] = -2d^{m 2 d - m 2 ) 2 sin 2 0. (28) 

Strictly speaking, in the spirit of the Ginzburg-Landau expansion it is not consistent to keep one C(M 4 ) term and 
drop all 0(M 2 ) terms. Guided by the high-density results we need to assume that the particular Ginzburg-Landau 
coefficient in front of the term (|28[) is large enough to compensate for the suppression due to the quark mass. 

We can now collect all terms in Eqs. (fT9|) . ([20)) . (j21[) . ([28l) and set m u ~ ~ to write the diquark potential as 



ax + a 2 m s cos <j> 2 /3x + /3 2 m s cos (/) - fi 2 sin 2 cj) 4 
\ld = ^ d "i 4 d ' 

where we have written the prcfactor of the term (|28[) as /i 2 , reminiscent of the effective kaon chemical potential ([2"T[) 
of the effective theory. 

D. Interaction potential 

For the interaction potential fl^d, the term of lowest order in the order parameters is 

Tr[d R dl§] + h.c. = 2d 2 [a. u + (a d + a s ) cos0] . (30) 



1 Alternatively, we may write the diquark field as 

/ 1 

d L = djj = di cos(<^/2) ie 1 ^ sin(</>/2) 

\0 ie-^ t sm(<f>/2) cos(<j>/2) 

including the effective kaon chemical potential fi. Then, the time derivative term, 

-Tridoid^doidndl)] = -2/i 2 d 4 sin 2 <f, , 

produces the structure of Eq. H28I I. This way of introducing a meson chemical potential in the time dependence of the condensate is 
done for instance in the context of meson condensation in nuclear matter W2H44I . 
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This term is anomalous since it is not invariant under U(1)a- Without kaon condensate it has already been considered 
in Ref. [Til We see that, in contrast to the pure diquark terms, the kaon condensate appears even in the absence 
of a strange quark mass. Including one more power in $ yields the following terms which are invariant under the full 
symmetry group G, 



Tr[di4 + d R d R }Tr[&<f>] = 6Tr[d L 4$$t] = 6Tr[d R d R $&] = 6d 2 {a 2 u + aj + a 2 s ) , 
dct^Tr^LG^* -1 ] + n - c - = 2d 2 [a d cr s + a u (cr d + cr s ) cos (j)] . 



(31a) 
(31b) 



We include mass terms up to 0(M ) which arise from replacing a chiral field $ in the above expressions (obviously, 
replacing a field dLd R in these terms does not yield interaction terms), 

Tr[d L d[ + d R d R ]Tr[M^] + h.c. = 6(Tr[d L d f L M&] + h - c -) = ^d 2 (m u a u + m d a d + m s a s ) , (32a) 
tabctijkM ai § b j(d L d) R ) c k + h.c. = 2d 2 {m d a s + m s a d + [m u {a d + a s ) + a u {m d + m s )] C0S(/>} . (32b) 

Estimates of the prefactor of the d 2 a 2 term suggest this contribution to be negligible [l9|, [2(|. We shall thus focus 
only on the d 2 a terms, as it was done for the three- flavor case in Refs. [ljl Hoj]. Then, collecting the terms in Eqs. 
pop . ()32|) . and again setting m u ~ m d ~ 0, <j u = a d = a s = a, the most general form of the interaction potential can 
be written as 



1 + 2 cos < 



+ m s (71 + 72 cos i 



d 2 a. 



(33) 



For 



we recover the interaction term — jd 2 a from Refs. [19l . |2C 



E. Full potential 



We can now put the contributions (fl8|) . (|29| . (|33j) together to obtain the full Ginzburg-Landau potential, including 
mass corrections and the meson condensate in the approximations discussed above, 

a\ + m s a 2 2 , c i + m sC2 3 b 4 
il = a m s (T H d H a + -a 

ol\ + a 2 m s cos 4> ^2 Pi + ferns cos 4> — fJ. 2 sin 2 4> ^4 



1 + 2 cos ( 



7" 



+ m s (71 + 72 cos i 



d 2 a. 



(34) 



In the remainder of this paper we discuss and evaluate this potential, i.e., we shall solve the stationarity equations 
for the three order parameters cr, d, and cj>, 



on on _ dfl 

da dd d(f) 



0, 



(35) 



in order to determine the ground state in the given parameter space. As written in Eq. Q34p there are 14 independent 
parameters. This is too unwieldy for a systematic study and therefore we shall work in several limit cases in the 
subsequent sections. In the next section we start with the case of a vanishing kaon condensate. 



III. MASS EFFECT WITHOUT MESON CONDENSATE 



For vanishing kaon condensate, = 0, all mass terms except for the linear a term are simply numerical corrections 
to the Ginzburg-Landau parameters. We thus absorb these mass terms into new overall prefactors. This reduces the 
number of independent parameters and does not change the results qualitatively (we do not attempt do determine the 
complete quantitative effect of the strange quark mass in the full parameter space). Then we can write the potential 
O as 

n(tr, d) = a a + \a 2 - C -a z + + |d 2 + ^d 4 - 7 d 2 a . (36) 
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The minus signs in front of the cubic term <r 3 and the interaction term d 2 a are conventions which imply c, 7 > [2C| . 
We also assume b, (3 > which ensures the boundedness of the free energy. The other parameters ao, a, a shall be 
varied without further constraints. For notational convenience we have absorbed the factor m s into the definition of 
ao. Without this linear term, the potential ([31)| has been discussed in Refs. [l9|, [2(|, where the phase structure in the 
(a, a)-plane has been given. We shall discuss the phase structure in the three-dimensional parameter space (a®, a, a) 
and present cuts through this space for fixed values of ao, such that we reproduce the known results for the special 
case ao = 0. We shall sec that several analytical arguments used in the massless case regarding phase transitions and 
critical points can be used repeatedly also for the more complicated cases. For the complete evaluation of the phase 
diagram, however, we need to employ numerical calculations. 
Without kaon condensate there are two stationarity equations, 



= -r— = ao + aa — co + bo' — 7a 
oa 

on , 

= — = ad + (3d 3 - 2jda . 
oa 



(37a) 
(37b) 



For the determination of the ground state it is useful to also consider the second derivatives. The Hessian matrix of 
Q(d, a) is 



H 



d 2 n 

da 2 dadd 

d 2 n d 2 Q 

V ddda ~blP 



a-2c<r + 3ba 2 -2-fd 

-2 7 d a + 3/3d 2 - 27CT 



(38) 



For a solution of Eqs. (|3"T1) to be a local minimum, the Hessian, evaluated at this solution, must have positive 
eigenvalues. If the potential is bounded from below - which is guaranteed by b, f3 > - the solution of Eqs. (|37|) with 
the lowest free energy yields the ground state, i.e., the global minimum, unambiguously. In other words, if we find 
a stationary point of the potential with a negative eigenvalue of its Hessian, then we must also have found another 
stationary point with lower free energy. Since the potential is bounded from below, the stationary point with lowest 
free energy must have a positive definite Hessian. Consequently, the second derivatives are strictly speaking not 
necessary for a stability check. We shall see below, however, that they yield a useful alternative method to determine 
the phase transition lines, some of which can be computed in an relatively simple way by using the eigenvalues of the 
Hessian. 

Employing the terminology of Refs. 0,H3]j we distinguish the following phases, 



NG phase : a + , d = , 
COE phase : a ^ , d ^ . 



(39a) 
(39b) 



In the Nambu-Goldstone (NG) phase there is a chiral condensate but no diquark condensate, while the coexistence 
(COE) phase is a color superconductor in the CFL phase with a nonvanishing chiral condensate. One might dis ting uish 

This was termed the normal (NOR) phase in Ref. 



19. 2C 



two more cases. First, both condensates may vanish. 
Clearly, in the presence of the linear a term, we cannot have two vanishing condensates, as we see from Eq. (|37ap . 
Only for a = this case occurs. It thus appears as a special case of the NG phase and is realized on the a = plane 
in our three-dimensional (ao, a, a) phase space. In all other cuts for fixed ao there can only be an approximate NOR 
phase with a small but nonzero value of a. A similar discussion applies to the phase with nonzero diquark condensate 
but vanishing chiral condensate, i.e., a "pure" CFL phase. With Eq. (|37ap and a = 0, the diquark condensate becomes 



d 2 = ^. (40) 

7 

Since we require d 2 > 0, a CFL phase is possible for ao > 0. Then, Eq. (|37b[) yields a constraint for the parameters, 



-a — . 

7 



(41) 



For a given ao, this is simply a straight line in the (a, a) plane where a vanishes. This line is of no particular interest 
and we obtain it as a limit case of the COE phase. Therefore we can restrict ourselves to the two phases given in 
Eq. (|39[) . First we shall discuss them separately in order to determine first-order phase transitions and critical points 
within the respective phase. Then we compare the free energies of the phases to obtain the phase diagram. 
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FIG. 1: Illustration of the appearance of a critical point for a free energy with a linear, quadratic, and positive quartic term 
in the order parameter, such as in Eqs. (|44p and ()51f) . The coefficients in front of the linear and the quadratic term change 
their sign across the two straight lines. For a negative quadratic term there is a first-order phase transition, indicated by the 
thick solid line. This line ends at the critical point at which both linear and quadratic terms vanish. Across the dashed lines 
there is no phase transition. This implies that one can smoothly connect two nontrivial minima of the potential by "going 
around" the first order line. In our context, these two phases are the NG phase and the (approximate) NOR phase or the 
COE phase and the (approximate) CFL phase. For these two cases, the location of the critical line and the critical point in 
the Ginzburg-Landau phase diagram are shown in the left and right panels of Fig. [2] 



A. NG phase 



In the NG phase, Eq. (|37b[) is automatically fulfilled, and Eq. (|37a[) yields a cubic equation for er, 

= a Q + aa - c<r 2 + ba 3 . (42) 

The general solution of this equation is lengthy and not very interesting. Throughout the paper, the determination 
of the phase structure requires solving cubic equations and insertion of these solutions into the free energy. This is 
best done on a computer and hence we do not attempt to derive all details of the phase structure analytically. Many 
features of the phase diagram, however, can be obtained in a simple way and understood on the basis of symmetry 
arguments as we explain in the following. 



By defining 



we can write the potential of the NG phase as 



t = *-±, (43) 



^ng(t) = Q a + a*T+^T 2 + -t\ (44) 



such that the cubic term has been eliminated. We have abbreviated 

„2 „4 



clqc ac 

na 55 ^ + 18^-108^' (45a) 



ac 2c 3 
3b ~ 27b 2 



r 2 



The cubic equation (|42j) then simplifies to 



a c = a - — . (45c) 
3b 



= a^+a c r + br 3 . (46) 

The phase structure within the NG phase is now determined by the signs of the linear and quadratic terms in the 
potential. Therefore, we need to determine the lines <2q = 0, a c — 0. Their intersection represents a critical point at 
which a first order phase transition line ends. This is explained for a general potential of the form (|44|) in Fig. [TJ In 
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FIG. 2: Left panel: first order phase transition line (thick solid line) and critical point in the (a, ao) plane within the NG phase, 
see Fig. [T]for a general explanation. The plot shows the coordinates of the critical point and of the intersections of the lines 
rag = 0, a c = with the coordinate axes, derived from Eq. 1)44[> . Since the NG phase is given by the chiral potential Q$ only, 
the coordinates depend on the Ginzburg-Landau coefficients b and c, and not on a, j3, and 7. In the three-dimensional (ao, a, a) 
diagram, the critical point thus becomes a critical line parallel to the a-axis where a first-order phase transition surface ends. 
Right panel: analogous scenario for the COE phase in the (a, a) plane (note different vertical axes of the two plots!), derived 
from the potential (|51|) . This critical point in the COE phase is present only due to instanton effects and vanishes from the 
phase diagram for 7 = 0, leaving behind a first-order line which does not end. 



the left panel of Fig. [5] we specify the location of the critical line and the critical point of the NG phase in the (a, ao) 
plane (the NG phase is obviously independent of a). 

In the case of the NG phase, the stability criterion from the Hessian ([55)) is very simple. The positivity of its 
eigenvalues is equivalent to the conditions 

a>2c<7- 3ba 2 , a > 2-/a . (47) 

The curve a = 2ja(a) in the (a, a) plane with a(a) being the solution to the cubic equation (|42p will turn out to 
coincide with the second-order phase transition lines between the NG and COE phases, as we demonstrate in Fig. [5] 
(The first condition, apart from being a crosscheck for the stability of the NG phase, is of no further interest to us.) 



B. COE phase 

For the COE phase we use Eq. (|37b[) to write the diquark condensate as a function of the chiral condensate, 

*(*) = ^p. (48) 

(In the numerical evaluation, one has to use this relation to check whether for a given solution a, the diquark 
condensate is real, d 2 > 0.) Inserting this into the potential ([55)1 yields 

ft COE [M»] =-^ + (ao + 7> + y<T 2 -^ 3 + ^ 4 , (49) 

where we have used the notation from Ref. [20l |. 

* _ "7 * _ 2 7 2 

7= T , a=a- — . (50) 

We can now proceed as for the NG phase. With the new variable t from Eq. we can write the potential as 

Ocoe(t) = ^ + ttc + (a + 7c>+f r 2 + h -r\ (51) 
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FIG. 3: Phase diagrams in the (a, a) plane without strange quark mass effects and without meson condensation, without (left 
panel) and with (right panel) anomalous effects. First-order phase transition lines are solid, second-order lines are dashed. For 
each phase we have indicated its global symmetries. They explain why the first-order phase transition line does not (left) and 
does (right) end at a critical point (see text and Fig. [4] for detailed explanations). The vertical (first-order) phase transition 
line separating NG from NOR (left and right) and COE from CFL (left) is located at a = 2c 2 /(9b). For the coordinates of the 
critical point (right) see right panel of Fig. [2] We have set the Ginzburg-Landau parameters to b = 1.2, c = 0.5, /3 = 1.6 for 
both plots, and 7 = (left), 7 = 0.1 (right). This particular choice is completely irrelevant for the topology of the left plot. 
For the right plot, one slightly different topology can be obtained in a different class of parameter values, where the phase 
transition between NG and COE is of first order along a certain piece of the transition line, see right panel of Fig. 2 in Ref. 
[20l ] . We translate the scenarios shown here into possible topologies of the QCD phase diagram in Fig. [6] 



where all ao's are written explicitly. A nonzero ao gives additional contributions to the constant and linear terms in 
t. All other terms are identical to the massless case [2(| (remember that we have absorbed the mass terms in the 
overall coefficients unless they have produced new structures in the order parameters), 





a 2 


* 

7 c 


a*c 2 
+ 18& 2 


4/3 




* 


= 7* + 


* 

a c 


2c 3 


7c 


~3b ~ 

c 2 


27b 2 ' 


«* 


= a* — 


3b ' 





c 4 



(52a) 



The stationarity equation is thus 



IO80 3 ' 

(52b) 
(52c) 

= a Q + 7* + a*r + br 3 . (53) 

Again wc can determine the transition line for the first-order transition within the COE phase and the corresponding 
critical point. We show the coordinates of this line in the right panel of Fig. [2l It is obvious that the critical point is 
a consequence of the anomaly: if the anomalous term vanishes, 7 = 0, the critical point disappears from the phase 
diagram because its a coordinate goes to —00 or +00, depending on the sign of ^fp- + ao (while its a coordinate 
remains finite). We also see that a negative ao shifts, for a given anomaly coefficient 7, the critical point towards 
larger values of a while keeping the a coordinate fixed. Remember that for both plots in Fig. [2] we only know that 
the respective phases are a stationary point of the potential, and we have not yet determined the global minimum. 
And indeed, for negative ao with sufficiently large modulus, the a coordinate of the critical point is shifted into region 
where the NG phase, not the COE phase, is the ground state, i.e., the critical point has disappeared. We demonstrate 
this scenario in the next subsection in Fig. [5j 



C. Ginzburg-Landau phase diagram, symmetries, and conjectured QCD phase diagram 

The next step is to compare the free energies of the NG and COE phases. We start with the case ao = in order to 
reproduce the results from Refs. [l9|, [2(|. We refer the reader to appendix A of Ref. [2(j for an analytical derivation 
of some of the phase transition lines. The results without and with anomalous term are shown in Fig. [3J 
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FIG. 4: (Color online) Diagram (left): topological space of U(1)l x U(1)r with generators aL,ctR £ [0, 2n]. Opposite sides of 
the square have to be identified (in particular, all four corners of the square correspond to the unit element). The resulting 
torus can also be parametrized by olr = {otR + c*l)/2 £ [0, 27r] and ola = (&r — q l)/2 £ [— 7T, 7r]. Table (right): symmetries of 
the free energy Q and the NG, COE, and CFL phases without and with axial anomaly. Thick (red) lines and points indicate 
group elements under which the free energy and the respective phases are invariant. For instance, in the first row, the points for 
the COE phase are obtained as the common subset of the points of the NG and CFL phases (and all sets of points are subsets 
of the points for Q). Only with anomaly, i.e., only after restricting the symmetry of Q from U(1)l x U(1)r to U(1)b X Za(G), 
the symmetries of COE and CFL are identical, allowing for a smooth crossover between these two phases. 



Since we have set the linear term in the chiral potential to zero, there is a NOR phase in both scenarios of the 
figure, where the chiral condensate vanishes exactly. If we switch off the anomaly, 7 = 0, the interaction term oc d 2 a 
vanishes. Since this is the only interaction term in our approximation, the chiral and diquark parts of the free energy 
decouple, resulting in a very simple phase diagram (left panel). In this phase diagram COE and CFL phases are 
separated by a first order line which does not end at a critical point (vertical solid line). This is a consequence of the 
global symmetries of the two phases. As indicated in the figure, the symmetries of COE and CFL are different without 
the anomaly. In other words, adding a chiral condensate to the CFL phase changes the (discrete) symmetries of the 
phase. This can be seen from the transformation properties of the order parameters in Eqs. (U]) and (J5]). The chiral 
condensate $ = diag(cr, a, a) and the CFL diquark condensate spontaneously break the chiral group SU (3)l x SU (3)# 
down to the vector subgroup of simultaneous left- and right-handed rotations SU(3)v- From Eq. (j4|) we see that any 
chiral condensate is also invariant under U(1)b transformations and under transformations of an axial subgroup 
Z^(2). However, this Z^(2) is contained in U(1)b- To see this it is helpful to consider the group J7(1)l x U(1)r as 
a topological space, in this case a torus, on which the discrete subgroups are sets of discrete points. We show this 
geometric picture in Fig. 2] which illustrates and facilitates the discussion of the discrete subgroups, in particular since 
we switch repeatedly between the bases of left- and right-handed rotations vs. axial and vector rotations, which can 
be confusing without this illustration. 

In contrast to the chiral condensate, the diquark condensates d^, da break U(1)b spontaneously. Hence there 
must be a true phase transition separating the NG and NOR phases from the COE and CFL phases. This statement 
is independent of the anomaly. The anomaly becomes important for the difference between COE and CFL phases. 
In the absence of instanton effects, the diquark condensates are invariant under independent sign flips of left- and 
right-handed quark fields, Z^(2) x Z#(2). This discrete group is broken by the chiral condensate; therefore, the COE 
phase, containing both order parameters, is only invariant under the common subgroup of U(1)b and Z^(2) x Zr(2). 
This is the group of simultaneous sign flips, Z(2). As a result, without axial anomaly COE and CFL have different 
residual symmetry groups, see first row of the table in Fig. [4] 

The anomaly reduces the axial symmetry U(1)a of the potential to Z^i(6), see Eq. (fT5)) . The group U(1)b X Z^(6) 
is represented in the first panel of the second row in Fig. 2] The anomaly does not affect the residual group of the NG 
phase. However, the residual group of the CFL phase is reduced since only Z(2) C Z^(2) x Z#(2), not the entire group 
Z/,(2) x Zfl(2), is a subgroup of U(\)b x Za(6), as can be seen geometrically in Fig. [4] Therefore, if instanton effects 
are taken into account, CFL is invariant only under the group of simultaneous sign flips Z(2). Now the addition of 
the chiral condensate does not further reduce this group, and the residual groups of COE and CFL become identical, 
see second row of Fig. 0] This allows for a smooth crossover between these two phases, and thus the first-order line 
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FIG. 5: (Color online) Ginzburg-Landau phase diagrams in the (a, a) plane with quark mass effect in the chiral potential 
(without meson condensation). Solid lines are first-order phase transitions, thick (black) dashed lines are second-order phase 
transitions. We have chosen various values of the linear coefficient cio of the chiral condensate and fixed all other parameters 
as in the right panel of Fig. [3] (in particular, 7 = 0.1). The parameter ao is given in units of ^frr because for ao < — ^ ne 
critical line which separates NG and approximate NOR phases vanishes and there is a crossover between these two phases. 
The thin (red) dashed line is the curve a — 2'y<j(a) with a(a) given by the solutions of the cubic equation (|42[) . This curve is 
obtained from the stability criterion of the NG phase and coincides with all second-order phase transition lines. 



between COE and (approximate) CFL can end at a critical point. 

This expectation from symmetry arguments is borne out in the Ginzburg-Landau phase diagram, see right panel 
of Fig. [3] This diagram will serve as a basis for our extensions in the following. 

One might ask why the straight line of the first order phase transition in the COE phase terminates at the intersection 
with the second-order line and then continues as a nontrivial, non-straight, curve. From the right panel of Fig. [2] one 
could have expected the straight first-order line to be present wherever the ao + 7* = line lives in the COE region. 
The reason for this behavior is the condition d 2 > 0. In terms of the schematic potentials of Fig. [TJ the first order line 
separates two nontrivial minima. Here, for positive a, the second minimum is forbidden since it would imply d 2 < 0. 
Consequently, upon crossing the line ao + 7* — the ground state remains in the same local minimum although, if 
one had ignored the condition d 2 > 0, there would have be a second local minimum with lower free energy. 

Let us now switch on the effect of the linear term in cr, induced by a nonzero strange quark mass. The chiral group 
SU(3)l x SU(3)r is now only approximate, which allows for a smooth crossover between the NG and (approximate) 
NOR phases, see left panel of Fig. [2] More precisely, within our ansatz M = diag(0, 0, m s ) we leave part of the chiral 
group intact, namely the SU(2)l x SU{2)r associated with the u and d quarks. On the other hand, our ansatz 
$ = diag(er, cr, cr) includes chiral condensates for the u and d quarks and thus one might think that $ does break the 
chiral symmetry in the u, d sector spontaneously. However, as one can see for instance from Eq. (|16p. once we have 
set cr u = ad = o~s it makes no difference (up to a numerical prefactor which is unimportant for our purpose) whether 
we set m u = rrid = or not. Therefore, although we consider the simplified situation of vanishing u and d quark 
masses, the NG and NOR phases are no longer distinguished by symmetry in the presence of a nonzero m s . As a 
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FIG. 6: (Color online) Conjectured translations of the Ginzburg-Landau phase diagrams to the QCD phase diagram in the 
(/is, T) plane [45|]. We have indicated the global symmetries which are broken in the various transitions. Left panel: zero quark 
masses, no instanton effects, corresponding to the left panel in Fig. [3] Middle panel: zero quark masses, nonzero instanton 
effects. In this case the CFL and COE phases are no longer distinguished by symmetry and thus allow for a smooth crossover. 
Whether the critical point is indeed present in the QCD phase diagram cannot be decided from the Ginzburg-Landau study; 
the T = axis may or may not intersect with the first order line in the COE phase in the right panel of Fig. [3] Here we show 
the case where it does not. Right panel: nonzero quark masses, nonzero instanton effects, see Fig.JS] Here we know from lattice 
QCD that the first-order phase transition between NG and (approximate) NOR does not reach the fin = axis. 



consequence, we expect the first-order line between NG and NOR in the uq = plane to become a first-order surface 
in the three-dimensional (ao, a, a) phase diagram which ends at a critical line at arj = — ^fp' - ^ n °th er words, if we plot 

phase diagrams in the (a, a) plane for fixed values of ao, this first order line should be absent for all ao < — ^fp - ^ ne 
numerical evaluation confirms this expectation, see Fig. [5] Besides the phase transition lines we have also plotted the 
curve a = 2ja(a) (red dashed line), where cr(a) is given by the real solution of Eq. (|42j) . This curve is obtained from 
the stability criterion of the NG phase, see Eq. (|47|) . The main purpose of this curve is to provide a semi-analytical 
form of the second-order phase transition lines. In Fig.[5]they have been obtained independently by directly comparing 
the free energies of the different phases. 

Besides the disappearance of the critical line between NG and approximate NOR phases (and deformations of the 
transition lines which do not change the topology of the phase diagram and thus are not very interesting for our 
purpose) there is one more topological change upon varying ao- Namely, we see that for sufficiently small values of 
a the first-order line within the COE phase disappears and the phase diagram in the (a, a) plane consists of a sole 
second-order transition separating NG and COE phases. More precisely, upon decreasing ao the critical point moves 
towards larger values of a while its a coordinate remains fixed. Interestingly, at the same value of ao where the first 
order line between NG and approximate NOR phases disappears, the critical point sits on the a = axis. We can 
see this analytically from Fig. [2] Then, upon further decreasing ao, it approaches the phase transition line between 
NG and COE phases and disappears for values of ao below some critical value for which we do not have an analytic 
expression. We have thus found two different ways to make the critical point disappear: switching off the anomaly 
removes the critical point but leaves the first-order critical line, while going to (possibly unphysically) small values of 
ao removes the critical point and the critical line. 

It is interesting to speculate how the Ginzburg-Landau results translate into the QCD phase diagram. For a precise 
translation our simple approach is of course not sufficient. Firstly, the Ginzburg-Landau potential is an expansion 
in the order parameters and thus cannot account for the complete potential; the expansion is expected to fail far- 
away from second-order phase transitions where fluctuations of the phase of the order parameter can become more 
important than fluctuations of the magnitude of the order parameter. We have also neglected gau ge fi eld fluctuations 
which in fact drive the transition from a color-superconducting to the normal phase first order |25l |46| . Secondly, 
even if we assume the Ginzburg-Landau approximation to be valid for all temperatures and densities of interest, we 
would need the dependence of the Ginzburg-Landau coefficients on the baryon chemical potential \ib and temperature 
T. This dependence is not known within full QCD, since the relevant regions of the phase diagram involve strong- 
coupling effects, and lattice calculations are inapplicable due to the sign problem at finite [Lb- We may still conjecture 
a translation from the Ginzburg-Landau phase diagrams to QCD, see Fig. [6j In this figure, the left and middle panel 
correspond to the situation without strange quark mass, and thus to the Ginzburg-Landau diagrams in Fig. [3] 

For the interpretation of our Ginzburg-Landau results with strange quark mass, it is helpful to think of the QCD 
(//b, T) plane to be a complicated surface in our (ao, a, a) parameter space. We know from lattice calculations that, at 
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an = 0, the transition from the chirally broken to the chirally (approximately) symmetric phase is a smooth crossover 
[471 ]. Therefore, the temperature axis must not intersect the critical surface between NG and (approximate) NOR 
phases. This critical surface may then manifest itself as a critical line between NG and approximate NOR phases 
which ends at a critical point, see right panel of Fig. [51 Another logical possibility is the absence of this line j48j . 
which would for instance be realized if the whole (hb,T) surface were located "behind" the ao = — ^fp- plane. 

For the critical point at low T and large hb our analysis with finite strange quark mass has opened up a third 
possibly topology. Without mass, the point may, although always present in the Ginzburg-Landau phase diagram, 
be either outside the {hb , T) plane (formally, one can think of the point being located at negative T) or within the 
(fiB ) T) plane. With mass effect we have seen that the critical line that ends at this critical point may be absent in 
the Ginzburg-Landau phase diagram. Thus, if the (/xb,T) surface is located at sufficiently negative ao, there is no 
first-order transition within the COE phase. 



IV. INCLUDING MESON CONDENSATION 



We can now extend the results from the previous section by allowing for a nonzero kaon condensate </>. In principle, 
this requires to consider several additional independent Ginzburg-Landau parameters as we can see from the full 
potential (|34p . In this potential, the parameters «2, /?2, ^ 2 , an d 72 become relevant when we allow for nontrivial 
values of 4>. All of these parameters correspond to mass terms and one might, as a first approximation, neglect these 
terms. However, then the only nontrivial structure involving is the unsupprcssed d 2 cr cos <fi term, and there would 
be no kaon condensation at all (in principle, there could be a condensate at the fixed value of (f> = 7r). Therefore, we 
have to keep the fj, 2 sin 2 4> term in order to match our potential to the high-density effective theory, see discussion in 
Sec. Ill C[ but for simplicity neglect the terms proportional to a2, 82, and 72. We have checked numerically that the 
inclusion of these terms can indeed make a difference to the topology of the phase diagrams with kaon condensate. 
We shall come back to this issue in the discussion at the end of Sec. [V] 

Within this approximation, the only additional parameter compared to the previous section is /x 2 , and our potential 
becomes 

a ? c o b 4 a , 2 8 — u 2 sin 2 . l + 2cos0 , 9 
Q(a,d,4>) = a <7 + -cr 2 - -cr 3 + -<t 4 + -d 2 + - — ^d 4 - 7 - ^d 2 a. (54) 

A comparison with the potential of the high-energy effective theory (|25[) shows that one can consider the interaction 
term d 2 a cos <fi as an effective, dynamical mass term for the kaon. The boundedness of the potential requires the d 4 
term to be positive, which yields an upper bound for /j 2 , 

\x 2 <8. (55) 

The stationarity equations (|35[) are 

= 
= 
= 

We distinguish the following phases 



The first two phases are the same as in the previous section. Note that <j) only appears in the potential when d 
is nonzero. This is clear since without diquark condensation there are no kaons to condense. Therefore, the NG 
phase does not depend on 0, and we have not specified its values in Eq. (|57ap . The value of 4> distinguishes between 
the phases with nonzero CFL order parameter, COE and COE- if . Again, as discussed for the case without meson 
condensate below Eq. ([3U)l . the NOR and "pure" CFL/CFL-A' phases are obtained as special cases from the NG and 
COE/COE-A phases. They only exist in a two-dimensional subspace of the three-dimensional (a ,a,a) parameter 
space and thus shall not be further discussed in the following. 



ao + aa — ca 



an 

~dd 

on ,„ . 
-— = d sin ( 



, c, 1 + 2 cos (j> , 2 
be - 7 -dr , 



= ad+(B - n 2 sin 2 4>)d 3 - 2j- 



2 7 M ,2 

— a d cos< 

3 2 



3 

1 + 2 cos ( 



-da , 



(56a) 
(56b) 
(56c) 



NG phase : a ^ , d = , 
COE phase : cr^O, d^O, = 0. 
COE- A phase: a ^ , d, ^ . 



(57a) 
(57b) 
(57c) 
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A. COE-A phase 



In order to compute the phase structure in the presence of a kaon condensate, we need to compute the free energies 
of the three phases (|57|) . The results for the NG and COE phases can be taken from the previous section. Thus we 
only have to discuss the COE- if phase. It is convenient to express <fi and d as functions of a. Solving Eq. (|56c[) for 
cos <f) and inserting the result into Eq. (|56b|) yields 

"••^"l = 30^) • (58a) 

*w - (58b) 

These expressions can be inserted into Eq. (|56a[) to obtain an equation for a. Equivalently, we can insert them into 
the potential (|54p and then minimize it with respect to a. The potential becomes 

2 CJ* h 

SW-xo [<r, d(a), 4>(<t)} = - 4(^3"" 2) + + 7> + Y ff2 - + i ff4 > ( 59 ) 



l£ = 777 137 ' < = a - V ( 7 ~ + 72 ) ' ( 6 °) 



where 

«7 . = _ 2 7V_L i 

3(/?-M 2 )' 0/4-0 9 + ^ 

We see that the potential of the COE- if phase has the same structure as the one of the COE phase (|49|) . with 
modified coefficients 7* a*. We may thus proceed analogously to determine the critical point within the COE-A 
phase. With the variable r = a — c/(36) from Eq. |43|) we obtain the potential 

n C OE-K°(T) = ^ + n c ,, + (ao + 1 l,)r + ^r 2 + ^r\ (61) 



with 



7^c a*c 2 



" 4(/3 - fj, 2 ) + 3b + 18& 2 108& 3 ' (62a) 



7c% - 7; + ^ - ^ , (62b) 
c 2 

< M = a M ~ 7T, ' ( 62c ) 

in complete analogy to Eq. (|51[) and resulting in the stationarity equation 

= «o + 7^ + <^ r + br 3 . (63) 

After solving this equation, one has to insert the real solution(s) back into Eqs. (|55|) to check whether d 2 (a) > and 
— 1 < cos 4>(<j) < 1. If one or both of these conditions are violated, the given solution for a has to be discarded. 

Since the cubic equation has the same structure as for the COE phase, we obtain an analogous first-order line as 
shown in the right panel of Fig. [5] From a,Q + 7* = a* „ = we can compute the location of the critical point. Its 
coordinates turn out to be 

acoE-KO = Tb + — {j-^ + 72J . «coe-ko = — [a Q + — 2 ). (64) 

We can compare this to the critical point in the COE phase whose coordinates are read off from the right panel of 
Fig.EJ 

c 2 2 7 2 ( c 3 \ 

acoE = ^ + — , acoE = --^ + ^j . (65) 

Before we turn to the phase diagram we can give a semi-analytical expression for the phase transition line between 
the COE and COE- A" phases. Since this phase transition will turn out to be of second order, wc consider nonzero, 
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but very small values of </>, i.e., we approach the second-order phase transition from within the COE- A phase. Hence 
we can divide Eq. (|56c[) by sin </> and then set cos <fi = 1 to get a simple relation between a and d. With the help of Eq. 
(I58bl) . we eliminate d from this relation and obtain the value of a at the phase boundary between COE and COE- A , 

Using the relation between a and r from Eq. (|4"5)l we insert this into the stationarity equation which then 
only contains the Ginzburg-Landau parameters. One obtains a cubic equation for a, which has a very lengthy, but 
analytical, solution a(a) for the second-order phase transition line between COE and COE- A' phases. 

Our next goal is to compute the phase diagram including meson condensation. It is useful to start with the following 
two questions. 

(i) What is the fate of the critical point in the (a, a) phase diagram in the presence of kaon condensation? We have 
seen above that a first-order line ending at a critical point is possible in the COE and in the COE- A' phase. 
More precisely, there is a "would-be" critical point with coordinates given in Eq. (pS)) which, if the COE phase is 
the ground state at this point, is a true critical point, and there is a "would-be" critical point with coordinates 
given in Eq. (|64|) which, if the COE- A' phase is the ground state at this point, is a true critical point. This leaves 
us with four logical possibilities. There may be no critical point at all when both "would-be" critical points are 
covered by the "wrong" phases; there might be one critical point, either in the COE or COE- A' phase; or both 
"would-be" critical points are realized if they are covered by the "right" phases. This classification is very useful 
since the information about the critical points determines, to a large extent, the topology of the entire phase 
diagram. 

(ii) Is there a region in the parameter space (here we mean all parameters except for a and a) for which the COE 
phase is completely replaced by the COE- A phase in the (a, a) plane? This question is interesting in view 
of the quark-hadron continuity. Recall that the existence of the anomaly-induced critical point opens up the 
possibility to go, at least at zero temperature, smoothly from the COE phase to the highest-density phase, the 
(approximate) CFL phase. If we now introduce a meson-condensed phase we might introduce an additional phase 
transition, separating the COE from the COE- A phase. This phase transition cannot end at a critical point 
because the kaon condensate breaks strangeness conservation U(l)s which is an exact symmetry of QCD, i.e., 
the CFL phases with and without kaon condensation have distinct residual symmetry groups. Only if we take 
into account the weak interaction which breaks flavor conservation, this line can end at a critical point. Another 
way of saying this is that through weak interactions the Goldstone mode associated with kaon condensation 
receives a small mass which has been estimated to be of the order of 50keV [49]. Here we do not consider 
such small E/(l),g-breaking terms and thus whenever both COE and COE- A phases are present in the phase 
diagram, they arc separated by a true phase transition. Such an additional phase transition could be avoided if 
the COE phase is completely replaced by the COE- A phase. 



B. Critical point (s) with meson condensate 



In the following we shall always use fixed values for the parameters 6, c, and /3 with respect to which the results are 
insensitive. (We do not have a rigorous proof for that but we have made several numerical checks which suggest that 
variation of 6, c, and /3 does not change our main conclusions.) To answer question (i) we thus determine for each 
value of the parameter set (ao,/i,7) the number of critical points in the (a, a) plane. To this end, we compute the 
ground state at the two "would-be" critical points (|64p and (|55|) and check whether the COE- A phase is the ground 
state at the point and whether the COE phase is the ground state at the point (p5|) . The result is presented in 
Fig. [7] in the (ao, /z) plane for two values of 7. 

The most obvious observation is that for large /i it is more likely to find the critical point in the kaon-condensed 
phase. The reason is simply that with increasing fi the kaon-condensed phase covers more and more space in the 
phase diagram which supports our interpretation of \i as an effective chemical potential. For sufficiently small values 
of ao there is no critical point at all (except for very large values of fi, just below its maximum value, for which there 
is a critical point in the COE- A phase). The reason is the same as already discussed for the COE phase in the lower 
right panel of Fig. [5] both points (|64p and (|55|) are in a region where the NG phase is the ground state, and thus they 
are not realized. Fig. [7] shows that the region without critical points is shifted to smaller values of ao for increasing 
instanton effects, parameterized by 7. We also see that with decreasing instanton effect, it becomes more likely to 
find the critical point in the meson-condensed phase. 

The regions with a critical point in the COE phase are not separated from the regions with a critical point in 
the COE- A phase by a one-dimensional line. We rather find a two-dimensional region in the (ao,/z) space where 
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FIG. 7: (Color online) Classification of parameter regions in the (ao,(i) plane for two different values of 7 according to the 
number of critical points (CP) in the (a, a) plane. Recall that ao is the coefficient of the linear term in the chiral condensate, 
induced by the strange quark mass, fi is the effective kaon chemical potential (bounded from above by fj, = /3 1 / 2 beyond which 
our Ginzburg-Landau potential becomes unstable), and 7 parametrizes the strength of instanton effects. The plots show four 
qualitatively different cases each of which is represented by a diagram in Fig. [8] corresponding to the four marked points 
in the right panel. In the dark grey regions there are, according to our numerical algorithm, two critical points. However, 
the corresponding phase diagrams show that one of them seems to lie on top of the phase transition line between COE and 
COE- if °, see upper right panel of Fig. [5] 



our numerical algorithm finds two critical points, one in each phase. However, a closer look reveals that in this 
region (dark grey in Fig. [7]) the point (|64[) seems to lie exactly on top of the phase transition line between COE and 
COE-A' . In other words, moving within the dark grey areas, the critical point seems to "drag" the second-order 
phase transition line which is attached to it. Fig. [7] also shows that this two-dimensional region is "squeezed" to a 
point in two instances. The first instance is at /i = (I/?) 1 / 2 , ao = — ^fp^ From Eqs. (|64|) and (|65|) we see that for 
these parameter values the two potential critical points coincide and sit at a = 0, i.e., on the a axis. On the other 
hand, Eq. (|66[) shows that for \i = (I/3) 1 / 2 the phase transition line between COE and COE- if is identical with the 
a axis because in this case a can only be finite if a = 0. Consequently, both points (|64]) and (|65|) coincide and lie 
on the phase separation line. Now, keeping y fixed and varying ao will keep both points together but move them 
away from the a axis, while the phase transition line remains unchanged. Hence, depending on which direction in ao 
one takes, a critical point will appear either in the COE (smaller ao) or the COE- if (larger ao) phase. The second 
instance is at positive ao, and the scenario is quite different here. Now the two points (|64|) and (|65|) do not coincide. 
Nevertheless, for the given parameters (for 7 = 0.1 we read off y, ~ 0.25^/ 2 , a ~ 0.85^ ) both points lie on the 
second-order phase transition line between COE and COE- if , and again varying the parameters by an arbitrarily 
small amount creates a critical point in one or the other phase. In contrast to the first instance this is a purely 
numerical observation. 



C. Phase diagrams with meson condensate 

The four qualitatively different scenarios obtained from Fig. [7] are represented in the four panels of Fig. El corre- 
sponding to the marked points in the right panel of Fig. [7] We have checked that our phase transition lines between 
COE and COE- A' in Fig. which are obtained by a brute-force comparison of the free energies, are reproduced by 
the curves a(a) discussed below Eq. (|66j) . Together with the thin dashed lines in Fig. [5] which reproduce the phase 
transition lines between the COE and NG phases, this means that we have a relatively simple, semi-analytical, form 
for all second-order lines. 

The first three panels (upper left & right, lower left) show the phase diagram for increasing y and all other 
parameters fixed. Not surprisingly, the COE- A phase covers more and more phase space with increasing ji. More 
interestingly, for all values of \x that are allowed in our approximation, a finite region of the COE phase without 
meson condensation survives. Here we have increased /1 up to 90% of its upper limit, but we have checked that this 
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FIG. 8: (Color online) Four Ginzburg-Landau phase diagrams including strange quark mass effects and kaon condensation, 
according to the potential (|54[) . The parameters b, c, /3 and 7 are chosen as in Fig. [5] while ao and fi are given in each panel 
and correspond to the marked points in the right panel of Fig. [7] The (red) arrows in the upper right panel indicate the a 
values for which the order parameters are plotted in Fig. [5] as functions of a. 



conclusion remains valid for all allowed values of [i. This seems to answer the above question (ii) with no, and meson 
condensation always appears to induce an additional phase transition line which does not end at a critical point. We 
have checked, however, that this statement depends on our approximation. Taking into account additional terms, for 
instance the ones proportional to oti, ,82, 72 in Eq. (j34p . it is possible to find regions in the parameter space where 
the COE- A' phase completely eliminates the COE phase from the phase diagram. 

The lower left panel shows that for large [i the COE phase survives in two disconnected regions. From the translation 
of the (a, a) plane into the QCD phase diagram, as discussed in Fig. [51 we can expect the T — axis to pass through 
the larger of these two regions, on the left-hand side of the first-order transition. The smaller strip on the right-hand 
side can be expected to be passed upon heating up the CFL phase, in agreement with NJL model calculations [H|. 
In both regions it is interesting to check whether a less symmetric color-superconducting phase than CFL becomes 
favorable. We discuss this possibility in the next section where we include the 2SC phase in our calculation. 

The upper right panel is representative for all points in the dark grey area in Fig. [7] Although our numerical 
algorithm shows that the COE- A phase is the ground state at the point ([M]) . this appears not to be a critical point 
since there is no first-order line attached to it. Numerically we find that the second-order phase transition between 
COE and COE- A' in the vicinity of this point is very strong, i.e., the kaon condensate <f> develops a sizable nonzero 
value on a much smaller parameter region than it does further away from this point. This can be seen in Fig. [S] where 
we plot the order parameters a, d, and as a function of a for three fixed values of a, as indicated in the upper right 
panel of Fig. [8] The middle panel of Fig. [9] shows the behavior close to the point (J64J) . We see that also the curves 
for a and d are not smooth around this point. 

The curves for the order parameters also illustrate the first-order transition at large, but still negative, a (right 
panel) and its smooth version at small a (left panel). Translated to the QCD phase diagram, we can think of the 
latter, if present at all, as being closer to zero temperature as the former. In the case of the first-order transition, here 
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FIG. 9: (Color online) Order parameters a, d, and cj> for three values of a as functions of a, corresponding to the phase diagram 
in the upper right panel of Fig. [S] As a normalization for the chiral and diquark condensates we have chosen their maximal 
values in the selected section of the (a, a) phase diagram. These values <r max and d max are assumed in the lower left corner, 
i.e., at (a, a) = (0, —0.13). The meson condensate is normalized by its maximum value n. Left (right) panel: onset of kaon 
condensation and crossover (first-order phase transition) to a phase where the approximate Z_l(2) x Z_r(2) symmetry is restored 
happen at different values for a. Middle panel: these two transitions happen simultaneously. 



taking place in the COE phase, both a and d are affected significantly. After the transition, the chiral condensate 
goes to zero for large a. This is as expected because the phase at large a corresponds to the (approximate) CFL-K 
phase. For the crossover, here taking place in the CFL-K phase, we see that the diquark and meson condensates 
are not much affected, only the chiral condensate decreases smoothly but drastically. The location of this crossover 
is given by the continuation of the critical line, see Fig. [TJ 

We can rephrase the main conclusion from Fig. [9] in the following concise way. There are basically two transitions: 
in the first, the chiral condensate goes to approximately zero; this is either a first-order transition or a crossover 
since the symmetry which gets restored is only approximate in the presence of the axial anomaly. The second is 
the onset of kaon condensation which is always of second order since the broken symmetry is exact (neglecting weak 
interactions). The transitions appear at two separate points in the left as well as in the right panel (and in different 
orders, comparing left with right). In the middle panel, they appear approximately at the same value of a, which 
seems to be the reason for the interesting, non-smooth, behavior in this case. 



V. UNDER WHICH CONDITIONS DOES THE 2SC PHASE APPEAR? 



So far our choice of the color-superconducting phases was inspired by high-density arguments. We have considered 
the CFL phase, which is present at asymptotically large densities, and the kaon-condensed CFL phase, which is the 
first adjustment of the CFL phase to the effect of a small strange quark mass within a weak-coupling approach. Our 
calculation, however, intends to shed light on the phase structure at moderate densities where less symmetric phases 
may appear, as discussed in the introduction. In this section we take into account one of these phases, namely the 
2SC phase. In the 2SC phase, all strange quarks as well as all quarks of one color, say blue, remain unpaired, i.e., 
Cooper pairs are made of red up/grecn down and green up/red down quarks. At weak coupling and parametricalh 
small strange quark mass the 2SC phase has larger free energy than either CFL or unpaired quark matter 
Phenomenological models such as the NJL model suggest that this may no longer be true at large coupling (Fll . 
and the 2SC phase (or variants thereof) may cover a region in the phase diagram between the low-density chirally 
broken phase and CFL. The appearance of the 2SC phase would clearly interrupt a possible quark-hadron continuity 
since 2SC does not break chiral symmetry and thus true phase transitions would be unavoidable between hadronic 
matter and 2SC and between 2SC and CFL. Building on NJL model calculations with t/(l)yi-breaking terms (27l.[53|. 
it has been argued that the 2SC phase indeed covers the potential anomaly-induced critical point for a wide region 
in the NJL parameter space [HI 2 . 

To get an idea about the possibility of a 2SC phase in our general Ginzburg-Landau formalism, we discuss the 
2SC phase in the simplest possible way. We shall not attempt to study the whole phase space with 2SC and meson- 



2 After completing this work, we have learned that in a new version (v2) of Ref. [28| an appendix containing a discussion of the 2SC phase 
in a Ginzburg-Landau approach, with some overlap to this section of our work, has been added. 
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condensed CFL. This would require the use of several additional Ginzburg-Landau parameters. Without meson 
condensation we shall be able, however, to make some general statements about the phase diagram including 2SC. 
On a qualitative level, we give some arguments about the addition of meson condensation at the end of this section. 
We also do not attempt to account for electric and color neutrality and beta equilibrium. These conditions are crucial 
for the appearance of non-CFL color superconductors since they impose constraints on the Fermi momenta of the 
various quark species. So far our phase diagrams have included only the CFL and CFL-if phases, and are thus not 
expected to be affected much by these constraints because the symmetric CFL pairing pattern ensures, at least at 
T = 0, that the number density of all quark species is the same. For the phase diagrams including the 2SC phase, to 
be discussed in this section, the neutrality constraint is more important. In our general Ginzburg-Landau approach, 
however, electric (or any other) charge cannot be defined. This is only possible if an explicit dependence of the 



Ginzburg-Landau parameters on the chemical potentials is assumed [21( , for instance using results from perturbative 
QCD or a phcnomenological model. Here we keep the parameters general and determine their range where the 2SC 
phase appears in the phase diagram. Neutrality and beta equilibrium can then be expected to yield relations between 
the parameters and thus define an "allowed" (neutral and beta-equilibrated) subspace of the full parameter space. 

Without meson condensate, the CFL order parameter is simply cLl = d R = diag(d, d, d), which is obtained from the 
more general order parameter © by setting <f> — 0. For the 2SC phase, the order parameter is d^ = d R = diag(0, 0, d) 
which describes pairing of only up and down quarks of two colors. In order to compare the free energies of 2SC 
and CFL we need to go back to the general Ginzburg-Landau terms. We shall for simplicity keep our assumption 
<I> = diag(cr, a, a) in both 2SC and CFL. As an alternative ansatz, accounting for the broken flavor symmetry, one 
might use the ansatz $ = diag(cr, cr, 0) for the 2SC phase. In this case, the potential becomes trivial because the 
d 2 a term that couples chiral and diquark condensates vanishes, and we have checked numerically that the 2SC phase 
appears nowhere in the phase diagram. For a more complete study one would have to include the d 2 <j 2 interaction 
term and/or allow for independent chiral condensates $ = diag(er u , er^, er s ). Here we proceed with the symmetric 
ansatz for $ and show that the 2SC phase appears in certain regions of the parameter space in accordance with 
physical expectations and with NJL studies. 

Within this ansatz, the chiral part of the potential 0$ can be taken directly from Sec. Inland is the same for CFL 
and 2SC. With the help of Eqs. (TlT)]) , (|2T)]) . and (f2"Tj) . we write the general form of the diquark part as 



Q d = ai [Tr[d L dl] + Tr[d R d R ]j + a 2 [Tr[d[Md R } + h, 

+ ft {(Tr[di4]) 2 + (T r [d R d R }) 2 } + p 2 Tr[d L dl}Tr[d R d R ] + ft {Tr[(d L d{) 2 } + T r [(d R d R ) 2 }} 

+ /3 4 Ti-[d R d[d L d R ] + ft (Tr[d f L Md R d R d R ] + h.c.) + ft {Tt[d R Md L d\d L ] + h.c.) 

+ ft (Tr[d{Md R }Tr[d R d R ] + h.c.) + ft (Tr[dlMd R }Tr[d{d L } + h.c.) , (67) 

where we have assumed the coefficients in front of terms that are related by an exchange of L and R to be identical. 
Here, a\ and a 2 are the coefficients for the d 2 terms without and with mass insertion, and ft, . . . , ft and ft, . . . , ft 
are the coefficients for the d 4 terms without and with mass insertions. For the d 2 a interaction terms (we neglect again 
the d 2 a 2 terms) we obtain from Eqs. ([50)) and 



n id = 7i (Tr[dfl4 $ ] + h - c -) + 72 (Tr[di4 + d R d R }Tr[M^} + h.c.) + 73 (Tr[d L d{M^} + h, 

M ai $ bj (d L d R ) ck + h.c.\ , (68) 

with 71 and 72, ■ • ■ ,74 being the coefficients for the terms without and with mass insertions. After performing the 
traces we can write the 2SC and CFL potentials as 

^cfl = fi* + (Qi+Q'2m s )d 2 + (ft CFL +ft CFL m s )d 4 + (7 1 +7 2 CFL m s )d 2 a, (69a) 
«2SC = n* + ( y + a 2 m s ) d 2 + (ft 2SC + pf G m s )d A + (y + 7 2 2S V.) d 2 a , (69b) 

with given in Eq. (fT8|) . We see that for the d 2 and d 2 a terms we know the ratio of the coefficients between 2SC 
and CFL for the mass-independent terms. Without mass corrections we can thus simply use 1/3 of the coefficient of 
the CFL phase to obtain the corresponding term for the 2SC phase. Including a small mass term then corresponds 
to a small correction to this 1/3. (In this spirit, it is irrelevant that we also know the ratio between the coefficients of 
the m s d 2 terms.) Such a statement is not possible for the d 4 terms where in general we need independent parameters 
even for the mass-independent terms. Therefore, the overall coefficient in front of the d 4 term in the 2SC phase must 
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FIG. 10: Ginzburg-Landau phase diagrams including the 2SC phase for two different values of ao (one should think of increasing 
the strange quark mass from left to right) and all other parameters fixed as in Fig. [5] The 2SC phase occurs after including mass 
corrections according to the potentials (|70[1 . for these plots we have chosen —am s = 7m s = 0.05. As a function of increasing 
| oto | (ao < 0) the 2SC phase first appears on the left-hand side of the phase transition between NG and COE, then additionally 
in a disconnected region on the right-hand side (left panel), before the two regions merge for sufficiently large |<zo | (right panel). 
COE(CFL) denotes the phase with nonzero chiral condensate and diquark condensate in the CFL phase (denoted simply COE 
in all previous plots), COE(2SC) denotes the phase with nonzero chiral condensate and diquark condensate in the 2SC phase. 
Had we not taken into account the 2SC phase, the (second-order) transition between COE and NG would have been between 
the shown COE(CFL)/COE(2SC) and COE(2SC)/NG transitions, see upper and lower right panels of Fig. [5] In this sense, 
the 2SC phase extends the color-superconducting area. 



be taken as a new parameter and cannot be expressed in terms of a single coefficient of the CFL phase. As a result, 
we can write the free energies in the following convenient way, 

£}cfi>, d) = a a a + V - ^ct 3 + h -a 4 + |d 2 + ^d 4 - 1 d 2 a , (70a) 

^2sc(^, d) = a cr + V - ^ct 3 + ^a 4 + | Q + am s j d 2 + ^d 4 - 7 Q + 7m s j d 2 a , (70b) 

where a = a\ + a2m s , ft = fti FL + /3^ fl to s , 7 = 71 + 7^ FL m s in order to reproduce the potential from Eq. ([36]) . The 
new coefficients a, ft, 7 depend on the parameters in Eq. (|69[) in an obvious, but irrelevant, way. 

We can now proceed analogously to the previous sections and determine the ground state of the system. Now we 
need to compare the NG, COE(CFL), and COE(2SC) phases, where COE(CFL) is the phase with coexisting chiral 
condensate a and diquark condensate in the CFL pattern (this phase was simply termed COE in the previous sections) , 
and COE(2SC) is the phase where a coexists with the diquark condensate in the 2SC pattern. The main question 
is under which conditions and where in the phase diagram the 2SC phase appears. The results of our numerical 
evaluation arc summarized in the following observations and in Fig. 1101 

• Our conclusions are insensitive to the value of ft. Note that ft is the only new coefficient in the 2SC potential 
that enters to leading order. It could thus have potentially complicated our results. We have checked, however, 
that variation of ft does not change any of the following statements and the topology of the phase diagrams 
shown in Fig. 1101 For this figure, we have thus simply chosen ft = ft. 

• If we neglect the mass corrections to the d 2 and d 2 a terms, i.e., a = 7 = 0, there is no 2SC phase in our (a, a) 
phase diagram for all values of ao. In other words, the mass effect through the linear term in the chiral potential 
oc ao is not sufficient to trigger the 2SC phase. More mass corrections are needed. 

• As soon as mass corrections a < or 7 > or both are switched on, there is at least one region in the (a, a) 
phase diagram for all ao where the 2SC phase is the ground state. (It is a numerical observation that only the 
given signs of a, 7 yield the results shown in Fig. [TO] different signs, i.e., a > and/or 7 < require sufficiently 
large positive values of ao for the 2SC phase to appear.) The 2SC phase appears in the expected regions of the 
phase diagram, separating the NG phase from the CFL phase. This is shown in Fig. [TO] for two different values 
of ap < 0. Increasing |ao| for negative ao increases the area which is covered by the 2SC phase. 
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• The phase transition from COE(CFL) to COE(2SC) is of first order. This is clear since two of the gap param- 
eters must change discontinuously in the transition from d^ = d,R = diag(d, d, d) to di, = dn = diag(0, 0, d). 
Additionally, our numerical results show that the chiral condensate is discontinuous at this transition. 

It is an interesting question whether meson condensation can prevent the 2SC phase from appearing. Comparing 
Figs. 151 and [TOl the answer seems to be no because the COE-if° phase does not reach areas close to the transition to 
the NG phase, and this is exactly where the 2SC phase lives. However, we have to remember that in our discussion 
of kaon condensation we have neglected all mass terms except for ao and the one associated with the kaon chemical 
potential. The omitted terms are exactly the ones that are needed to obtain the 2SC phase. Hence we would have to 
redo our analysis, including 2SC and kaon condensation and taking into account the mass corrections for the d 2 and 
d a terms in the potential (|34|) . This would introduce more independent parameters and without further constraints 
a systematic study would be very unwieldy. Therefore, we have only done some numerical calculations with selected 
parameters. These calculations show that if we take for instance one of the phase diagrams in Fig.[T0l there is a range 
of parameters for the mass corrections where the kaon-condensed phase does expel the 2SC phase. However, the 2SC 
phase is only expelled completely from the (a, a) phase diagram if the mass corrections oc a.2, P%, 71, 72 in Eq. (|34[) 
are of the order of or larger than the 0(mP s ) terms. In this case, our Ginzburg-Landau expansion becomes unreliable 
since we have neglected mass terms of higher order (except for the term cx /1 2 which is suggested to be relevant from 
high-density arguments, as explained). If we keep the mass corrections much smaller than the 0{m s ) terms, the 2SC 
phase, if it is preferred over the CFL phase in some region of the phase diagram, also appears to be favored (in a 
smaller region) over the CFL-if phase. 

VI. SUMMARY AND OUTLOOK 

We have studied phases of dense matter in a Ginzburg-Landau approach. Previous Ginzburg-Landau studies have 
shown that the axial anomaly may induce a high-density critical point in the QCD phase diagram, possibly leading 
to a smooth crossover between hadronic matter and color-flavor locked quark matter. We have explained in detail 
that the existence of this critical point is a consequence of the (discrete) symmetry of the CFL phase which - in the 
presence of the axial anomaly - is not changed by adding a chiral condensate. Our main goal has been to extend 
the previous studies by including a strange quark mass. We have discussed several different, although related, mass 
effects. 

Firstly, the strange quark mass introduces a term linear in the chiral condensate cr, say ao<r, which allows for a 
smooth crossover between the phases of broken and (approximately) restored chiral symmetry. This effect is most 
relevant for the high-temperature, low-density phase of QCD where there is indeed such a crossover between the 
hadronic phase and the Quark-Gluon Plasma, as we know from lattice calculations. We have shown that the term 
aoc is also relevant for the high-density critical point. The reason is the anomalous interaction term that couples 
the chiral to the diquark condensate d, say r yd 2 a. For sufficiently large values of |ao| (ao < 0) the first-order phase 
transition line which, for nonzero 7, ends at the high-density critical point disappears. As a result the transition 
between the ordinary chirally broken phase and the CFL phase is smooth everywhere. 

Secondly, a nonzero strange quark mass is expected to induce less symmetric color-superconducting phases. In 
high-density calculations, the CFL-K phase is the first phase that appears after going down in density from the 
asymptotically dense CFL region. We have introduced a kaon condensate as a relative rotation of left- and right- 
handed diquark condensates and have adjusted the Ginzburg-Landau potential to match the essential terms of the 
high-density effective theory. We have identified the region in the parameter space where the critical point has moved 
from the CFL into the CFL-A' phase and have determined the location of both possible critical points in the presence 
of a strange quark mass. In addition to a shift of the critical point, the kaon condensate introduces a true phase 
transition because it breaks strangeness conservation spontaneously which is an exact symmetry in QCD. 

Thirdly, we have discussed a more radical reaction of the system to a nonzero strange quark mass, namely the 
appearance of the 2SC phase. In previous studies in the Ginzburg-Landau approach, it has been shown that for 
infinitely large strange quark mass there cannot be a high-density critical point. The reason is that the 2SC phase 
does not break chiral symmetry and thus there must be a true phase transition between the 2SC phase with and 
without coexisting chiral condensate. Since we have included a nonzero, but finite, strange quark mass, we could 
study the competition between the 2SC and CFL phases under the influence of a nonzero chiral condensate. We have 
shown that the mass term a^a is not sufficient to favor the 2SC phase in any part of the phase diagram. Additional 
mass terms, which we have neglected in our discussion of the meson condensate, are necessary for the 2SC phase to 
appear between unpaired quark matter and the CFL phase. In a conjectured translation to the QCD phase diagram 
it seems that the 2SC phase appears "first" (i.e., for the smallest values of these mass terms) at low temperature. As 
a consequence, the smooth crossover at zero temperature between hadronic and quark matter would be disrupted by 
true phase transitions. The appearance of the 2SC phase is in agreement with recent NJL model calculations. 
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There are several possible extensions of our work. We have studied the competition between 2SC and CFL system- 
atically, but have only briefly discussed the competition between 2SC and CFL- if . For a more complete analysis it 
would be helpful to first find some constraints for the additional Ginzburg-Landau parameters. The potential prolif- 
eration of parameters has also led us to a simplified, flavor-symmetric ansatz for the chiral condensate. One should 
check in further studies how our results change with a more realistic ansatz. It would also be interesting to consider 
different meson condensates and possibly their coexistence. This is important since we do not know the masses of 
the CFL mesons at intermediate densities, and it may well be a different meson than the kaon which condenses in 
this regime. Furthermore, one should also take into account the requirement of electric and color neutrality. This 
was not an issue in our calculations with the CFL phase since this phase is automatically neutral (and a neutral kaon 
condensate does not change this). In the 2SC phase, however, the numbers of up, down, and strange quarks are not 
identical, and the existence and details of this phase (as for any non-CFL color superconductor) depend strongly on 
the neutrality constraint. 

More generally speaking, the model-independent Ginzburg-Landau approach, including possible extensions in the 
future, is helpful to gain insight into the QCD phase diagram at low temperature and large, but not asymptotically 
large, densities. Like NJL model calculations, however, it is far from being conclusive for the actual, full QCD situation. 



Therefore, it is important to also pursue other approaches such as improvements of perturbative calculations [54j or 
studies of dense matter in the astrophysical context and comparing properties of phases of dense (quark) matter with 
data from compact stars (5^4131 • 
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